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Abstract 

This paper is a short summary of the main results in the thesis [T]. 
Based on the P2P paradigm we construct a stochastic model for a live media streaming con- 
tent delivery network. Starting from the behavior of the out degree process of each node we 
analyze the evolution of the resulting random graph. For the out degree process we provide an 
iterative formula for the distribution probabilities and calculate an explicit expression for the 
transient distribution. Moreover we discuss the steady state distribution and the problem of 
local disconnections. 



1 Introduction 

Multimedia applications have experienced explosive growth in the last decade and have become a 
pervasive application over the Internet. A video streaming system distributes a video stream syn- 
chronously to a large population of viewrers. Today, one-to-many commercial streaming solutions 
are offered via content delivery networks (CDNs) such as Akamai, Limelight or VitalStream. Such 
content delivery networks are based on P2P {peer to peer) paradigm, i.e. all users (peers) contribute 
on distribution service. 

In Section 2 we present a model for a live media streaming content delivery network. The construc- 
tion of the model is based on P2P paradigm according on P],[3]. We consider that the number 
of nodes which lie in the network is managed by a stochastic process and give a stochastic rule to 
create connections among the nodes. We focus our analysis on the study of the evolution of the 
structure starting from the local behavior of in and out-degree of a single node (user). 
In Section 3 we discuss the local characteristics of the network in the case of zero mortality. In this 
case the process associated to the out-degree is a pure birth process. After having provided the 
distribution probabilities of this process and its subordinated chain, we show a relationship with a 
Bernoullian scheme for which the successes probabilities are dependent by the number of previous 
successes. This approach is very useful for recovering some interesting properties like mean, vari- 
ance and for recovering concentration and deviation bounds. 

In Section 4 we turn to the general case. In this case the evolution of the out-degree process of the 
generic node is equivalent to a problem in queuing theory, this process is known in literature as dis- 
couragement queue. Although from the theoretical point of view the problem of the determination 
of the transition functions has been solved [3], the explicit form of them is not present in literature. 
Assuming that the solution is represented by a Taylor series we obtain an exact formulation of the 
distribution probabilities being determined recursively, then we demonstrate the convergence of the 
series. Moreover for the related embedded chain we obtain an explicit expression of the transition 
functions. Furthermore we discuss the out-degree distribution of the generic node in the steady 
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state and analyze a mechanism for recovering the connections in the case of local disconnections. 
We propose a solution in terms of "extra job" for the root (master server) and give an estimate for 
the mean of this extra job. 

2 The Model 

We consider a Poisson process, with parameter A > 0, that counts the number of accesses to the 
network. Any access corresponds to one user who wants to benefit of the streaming service. We 
consider that the life-time of each user is exponentially distributed with parameter ^ > 0. Thus the 
number of nodes that lie in the network at some time t > is subordinated to an M/M/oo queue. 
Remembering that in the P2P context each node in the network participates in distributing the 
service to a subset of [future) nodes, we give the following rule for generating connections among 
the nodes: 

At time of the access each node sends a service request concurrently to all nodes already in the 
system. Any request will he accepted with probability where k is the out-degree of the node that 
receives the request. 

For each accepted request we generate a directed connection {directed arc) from server node to 
client node. 

The connection probabilities have been chosen considering that every node decreases its capacity 

[availability to service) while the number of its clients increases. 

When a node drops from the network we erase this node together all its arcs. 

We consider the first node as the root [master server) and suppose this node to have an infinite 
life-time. We observe that by construction the out-degree processes of the nodes are mutually 
independent and the number of in and out connections of a node is a stochastic process. 
The generated graph is directed and acyclic, in fact all arcs are directed so that no one of the nodes 
receives informations from the nodes which are arrived later. 

2.1 Recovering out-degree processes 

Let {N[t)}t>o be the process that counts the number of nodes, except the root, in the networks. 
We know that N[t) is, in particular, a birth and death process with parameters A; = A, Z > 0, and 
= 111, I > 1 and fiQ ~ 0. 

Let us consider a generic node that is in the network at time t > 0. From the local point of view 
also the out-degree process of the generic node is managed by a birth and death process. 
Since the birth rates depend only on the parameter of Poisson process and on the probability of 
acceptance of the request, we have A^ = j^, k > 0. 

As regards the mortality rates, let us observe that if the process N[t) has a backward transition 
the state of the system changes because one of the nodes will be erased together with all its arcs. 
This event corresponds to a backward transition of the out-degree of the nodes that were adjacent 
to the node which has been erased. 

Assuming that there are n nodes in the network at time t, for a generic out-degree process the 
mortality rate will be the product of the mortality rate of the process N[t), i.e. n^, and the 
proportion of the number of neighbors of the node in relation to the total number of nodes. Thus if 
the generic out-degree was in state k < n then /i/j = -^n/^ = fc/i. 

In conclusion being {X{t)}t>{) the generic out-degree process, let P[X{t) = k) = p[k,t) then for all 
k G Nq we have: 
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p'{k, t) = - i^-^—^ + k^lj p{k, t) + -p{k - 1, t) + i{fe>i}M(fc + i)p{k + 1, t) (1) 

with initial conditions p(fc,0) = 5kfl- 

In the fohowing we discuss separately the case = and ^ > 0. 

3 Case // = 

In this case the equation ([T|) is reduced to 

\ p'{k,t) = ^pik^l,t)~j^pik,t), if k^O 
with initial conditions p{k,0) = 5k,o. 

A direct use of the Laplace transform gives us the following result. 
Proposition 1. Let 

^(fc.*)=^EW+-^/(';')e-^ (3) 
■ j=i \ J / 

then Eq.(^ is the solution of Eq.(^ with initial condition p{k,0) = Sk.o- 

Many interesting properties of the network do not depend on the elapsed time between nodes arrival, 
i.e. the Poisson process which counts the number of nodes in the network at generic time t, but 
only on the number of the nodes. From this point of view if we consider the subordinated chain 
X : No No of the generic out-degree process with transition probabilities pfc^fc+i = = l—pk,k, 
k > 0, using generating function approach we have the following 

Theorem 2. Let {Xn}n>o the subordinated chain of the birth process {X{t)}t>o associated to 
Eq.^ with A = 1, if we put P{X„ = k) = pn,k then 



Pn+l.k = -kTlPn,k + \Pn.k~l, Vn>0, l<fc<n+l 

Pij =0, i > 0, J > i + 1 

P0,0 = 1, Pnfi =0, n > 1 

Furthermore 

2—1 ^ ' ^ ' 

As a consequence of Eq. Q we establish a useful Bernoullian scheme for the subordinated chain. 
Proposition 3. Let Z„ = X^'jLo ^i, n>0, where 

then Zn = Xn, where the probability distribution of Xn is given by Eq.^. 
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The above Bernoullian scheme can be interpreted as a special type of Polya urn: 



Consider an urn containing only one white hall and an arbitrary number of red halls. If we draw 
the white ball then we add a red hall, while if we draw a red ball we do not do anything. In both 
cases we reinsert the drawn ball and proceed to the next drawing. 

From the previous proposition and some resuhs due to MacDiarmid [S] we estabUsh the following 
deviation bounds. 

Proposition 4. Let {X„}„>o the process with probability distribution function as in Eq.^, it 
results that 



Furthermore 



P{Xn > (1 + e)E{Xn)) < cxp ( -^e^E{Xn) ] < e < 1; 

1 
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P{Xn < (1 - e)E{Xn)) < cxp ( ^le^E(Xn) ] < e < 1 



3.1 The in-degree process: mean and deviation bound 

Let us remember that when a node have accessed to the network it sends a request to the nodes 
already in the system to establish a connection. The accepted requests are subordinated to out- 
degree processes because any possible transition for a generic out-degree process depends on its 
associated urn drawing. 

Let us consider n + 1 — i trials, i < n, of the i-th urn associated with the corresponding out-degree 
process. 

If X~_^^ is the in-degree process of the generic (n + l)-th node then 

n 

Xn+1 = 51 ^"+1 
i=l 

We observe that this value will not change after the (n + l)-th time step, and remember that the 
event {i^^+i = 1} corresponds to have generated the arc i ^ n + 1. Moreover being the urn schemes 
mutually independent then '^^ Eq.® is a sum of independent Bernoullian r.v.'s. Now we can 

give an estimation of the mean of the generic in-degree process. 

Lemma 5. Let X^^^ as in Eq.^ then 

EiXn+i) - (6) 

Before giving a deviation bound for the in-degrec process we need to estimate the probabilities of 
success for the n + 1 — i, n > i, drawing of the z-th urn, i.e. of the event = 1}. 

Lemma 6. Let 0n_^i = P{Yjl_^_-^ ~ 1), with i = 1, ...,n then 

V2(n-i) + l 
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Using an approach based on self-bounding functions, see [5], we establish a two sides deviation 
bound for the generic in-degree process. 

Proposition 7. Let ^,7+1 Eq.(^ then the sequence ^,7+i {^,0)-self-bounding and we have 

FIX- - B(X-) > ,E{X-}) < oxp (- ''^'f'"' ) ~ cxp f-^j^j 

P(xr - Eixn < -.E(xr)) < .xp {-'-^) ~ «p (^^) 

3.2 Spanning trees 

Starting from the in-dcgrce process we establish a function that counts the number of spanning 
trees and recovers the mean. 

We know that the arcs i ^ i + 1, with i E N, exist with probability 1 (by construction) and there 
is at least a spanning tree for any generated graph. 

Let Gn+i be the graph at time step n + 1 and let T„+i be the number of spanning trees in Gn+i- 
Recalling that at time step rt + 1 we add the node n + 1 and, with their probabilities, the arcs 
i — )■ n + I, with i = l,...,n, the number of these generated arcs represents the in-degree of the 
n + 1-th node. 
Observing that 

n 
i=l 

where T^-^^^i — being „_)_i') the indicator function of the arc i — > n + 1, i.e. T^+i is 

the number of spanning trees which contain the arc i — > n + 1 , if it has been generated, then from 
Eq.® we have 

n j 

j=i 1=1 

Because for every fixed i the Y^, with n = i,i + 1, ... are negatively correlated Bernoulli r.v.'s and 
for every fixed n > 1 the y^, with i < n — 1, are mutually independent Bernoulli r.v.'s, from Eq.® 
computing the mean of the number of spanning trees in G„+i we have 

4 Case // > 

If ^ > the out-degree process of the generic node is a so called discouragement queue, see [7]. 
In the following theorem we provide an exact formulation of its distribution probabilities being 
determined recursively. 
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Theorem 8. Let 



{0(1) _ ciO) , . c.(0) 
g{k+i) ^ ^ ^^_^|fc) _ ^2gik-i) ^ > 0^ < fc < i + 1 
^fe = ITT + ^ ^ 0' " = v/f 
6eOT<? S"!"^^ =0 Vi > and ^^"^ =0 Vn > i, t/ien V/c > 



1=0 

is the solution of Eq.(^, with initial conditions 0) = (5fe,o- 

The above result is obtained making use of Taylor series and showing that the coefficient 5'|''^ are 
bounded by Bcssel numbers. For more details on Bessel numbers see [8]. 

As regards the probability distribution of the embedded chain of the process {X{t)}t>o associated 
to Eq.(IT]), proceeding by induction we establish the following result. 

Proposition 9. Let {pn,k)n,k>o the embedded process of the birth and death process {X{t)}t>o 
associated to Eq. ^^ it result that 



I Pn+i,k = T+feptiy^Pn.fc-i + if(t+2Kt+i')a^ P".fc+i' > l<k<n+l 

\ Pn,k = 0, Vfc > n, pofi = 1 

furthermore 

Pn,k = dkTf, ^ , n>0, 0<fc<n, n + k even 
Pn,k = 0, otherwise 

d'k = rii^O 

, Ti'"^ = EJ^=o '''^V/'^' ^+r'\ Vfc >Q,h> 1; ^ = 1, Vfc > 



4.1 Steady state 

Taking into account that a M/M/oo queue has an exponential rate of convergence to equilibrium, 
see [9] , for applications it is interesting to know the behavior of the degree distribution in the steady 
state. We know that in the steady state the number of the nodes follows the Poisson distribution 
with parameter ^, then we can recover the degree distribution of the generic node starting from 
the events = / > 0. 

Let us consider the steady state and label the nodes with the natural sequence choosing the index n 
for node position. From Eq.(IT]), conditioning to N = I, we have that the distribution in the steady 
state verifies 

= + ^k+iMk + 1) + (-Y^ - Mfc) (8) 

Moreover if fc > / — n > we have tt^ ; = because the out-degree of n-th node can not be greater 
than I — n conditioned to the event N ^ I. 



6 



Proposition 10. //tt^j satisfies Eq.(^ then for k < I — n and n < I with l,n = I, 2, 3, ... it results 




furthermore for the unconditional distribution tt^ it results 

, \ fc -A +00 ( X \l 

A\ e 1^ M-^ {p 

a) (fc!)2 ^ /I V'"" (^)m ^ 1 _ 

4.2 Recover connections from the root 

In the general case it is possible for a node to have in-degree at a generic time t. This generates 
a temporary inefhciency of the streaming service. In this case we create with probability 1 a 
connection from the root to this node and restore the service. We will not modify the out-degree 
process for the root, we consider the act of connection recovery only as an extra work for the root. 
Let us consider the steady state. We know that the number N of the nodes {except the root) 
is distributed in accordance with Poisson distribution with parameter ^. Moreover the number of 
connections that each node established to benefit of the stream service is related to transitions of the 
out-degrcc processes of the nodes that was early in the network. Because the out-dcgrce processes 
arc independent these numbers of (in) connections is a sum of BcrnouUian random variables. 
Let X~ the random number of in-arcs of the generic node, from the fact that all nodes have the 
same distribution for the life time and from the minimum order statistic, we deduce that the extra 
work for the root is: 

^ 1 

Wroot = E 7^ (9) 

i=l i 

By virtue of Anscombc transform and of a generalized HoefFding inequality [10] we can established 
the following result. 

Proposition 11. Let Wroot o,s in Eq.^ then 



E{Wroot) ~ 
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